Accepted for Publication in the Astrophysical Journal 

Preprint typeset using I4TgX style emulateapj v. 1 1/10/09 



ASSESSMENT OF STELLAR STRATIFICATION IN THREE YOUNG STAR CLUSTERS IN THE LARGE MAGELLANIC 

CLOUD. 

DiMITRIOS A. GOULIERMIS 
Max Planck Institute for Astronomy, Konigstuhl 17, 69117 Heidelberg, Germany 

DOUGAL MACKEY 
Institute for Astronomy, University of Edinburgh, Royal Observatory, 
Blackford Hill, Edinburgh, EH9 3HJ, UK 
Argelander-Institut fUr Astronomic, Rheinische Friedrich-Wilhelms-Universitat Bonn, 
Auf dem HUgel 71, 53121 Bonn, Germany 

YU XlN AND AND 
BOYKE ROCHAU 

Max Planck Institute for Astronomy, Konigstuhl 17, 691 17 Heidelberg, Germany 
Accepted for Publication in the Astrophysical Journal 

ABSTRACT 

We present a comprehensive study of stellar stratification in young star clusters in the Large Magellanic Cloud 
(LMC). We apply our recently developed effective radius method for the assessment of stellar stratification on 
imaging data obtained with the Advanced Camera for Surveys of three young LMC clusters to characterize 
the phenomenon and develop a comparative scheme for its assessment in such clusters. The clusters of our 
sample, NGC 1983, NGC 2002 and NGC 2010, are selected on the basis of their youthfulness, and their 
variety in appearance, structure, stellar content, and surrounding stellar ambient. Our photometry is complete 
for magnitudes down to mg^ ~ 23 mag, allowing the calculation of the structural parameters of the clusters, 
the estimation of their ages and the determination of their stellar content. Our study shows that each cluster 
in our sample demonstrates stellar stratification in a quite different manner and at different degree from the 
others. Specifically, NGC 1983 shows to be partially segregated with the effective radius increasing with 
fainter magnitudes only for the faintest stars of the cluster. Our method on NGC 2002 provides evidence of 
strong stellar stratification for both bright and faint stars; the cluster demonstrates the phenomenon with the 
highest degree in the sample. Finally, NGC 2010 is not segregated, as its bright stellar content is not centrally 
concentrated, the relation of effective radius to magnitude for stars of intermediate brightness is rather flat, 
and we find no evidence of stratification for its faintest stars. For the parameterization of the phenomenon of 
stellar stratification and its quantitative comparison among these clusters, we propose the slope derived from 
the change in the effective radius over the corresponding magnitude range as indicative parameter of the degree 
of stratification in the clusters. A positive value of this slope indicates mass segregation in the cluster, while a 
negative or zero value signifies the lack of the phenomenon. 

Subject headings: Magellanic Clouds - galaxies: star clusters - globular clusters: individual: NGC 1983, 
NGC 2002, NGC 2010 - Hertzsprung-Russell diagram - Methods: statistical - stellar dy- 
namics 



1. EVITRODUCTION 

Stellar stratification is a characteristic phenomenon of star 
clusters, which has been well documented for more than 50 
years. It is directly linked to the dynamical evolution of star 
clusters dChandrasekhar ll942l : ISpitzedl987t) . and specifically 
to the central concentration of the massive stars as the clus- 
ter relaxes dynamically, a phenom enon known as mass seg- 
regation (for detailed rev iews see Lightman & Shapiroll 19781 : 
iMevlan & Heggielll997h . This behavior in massive Galactic 
globular clusters (GGCs) is well understood and the observed 
mass segregation is explained as a result of their dynamical 
relaxation in a timescale, freiax, much longer than their ages, 
r, as derived from stellar evolutionary models. This dynami- 
cal mass segregation, thus, forces more massive stars to sink 
inwards to the center of a cluster through weak two-body in- 
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teractions. 

However, segregation of massive stars is observed also 
in young star clusters, such as the ^ 30 Myr-old clus- 
ter N GC 330 in the Small Mageflanic Cloud (iSirianni et alj 
|2002t), or even the - 1 Myr-old Galactic starburst NGC 3603 
(Stolte et al. 2006). Such systems, supposedly not sufficiently 
old to be dynamically relaxed, cannot be included in the "tra- 
ditional" picture of dynamical mass segregation as it occurs in 
GGCs. For these cases the phenomenon has been - rather ten- 
tatively - characterized as primordial mass segregation, based 
on the suggestion that this stellar stratification is due to the 
formation of the stars in, or near, the central part of the clus- 
ters rather than being dynamically driven 

Mass segregation at early stages of clustered star forma- 
tion is indeed predicted by theoretical studies, which sug- 
gest that the positions of massive stars for rich young clus- 
ters of r << fieiax cannot be the result of dynamical evolu- 
tion jBonnell & Davies|[T998h . Moreover, the appearance of 
gas in such clusters reduces further the efficiency of any dy- 
namically driven mass segregation. According to theory, mas- 
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sive protostars can be formed at the central part of the pro- 
tocluster through cohesive collisions and dissipative merging 
of cloudlets, which lead to extensive mass segregation and 
energy equipartition among them (Murray & Lin 1996), or 
through accretion from a distributed gaseous component in 
the cluster, which leads to mass segregation since stars lo- 
cated near the gravitational center of the cluster benefit from 
the attraction of the full p otential and accrete at higher rates 
than other stars (iBonnell & Bate 2006). 

Both of the aforementioned massive star formation mecha- 
nisms predict primordial mass segregation through early dy- 
namical processes at the central part of the protocluster. The 
significant difference of these processes from the subsequent 
energy equipartition of the stars in the formed cluster is that 
the latter, and consequently dynamical mass segregation, is 
independent of the initial cluster conditions, while the former, 
and therefore primordial mass segregation, is directly con- 
nected to these conditions. In order to achieve a complete 
understanding of the phenomenon of primordial mass segre- 
gation, many related studies are focused on clusters at very 
early stages of their formation, including a variety of types of 
stellar systems from loose stellar associations and open clus- 
ters to compact massive starbursts. 

The Magellanic Clouds (MCs), being extremely rich in 
young star clusters (e.g., iBica & Schmitj 119951 : iBica et all 
119991) . offer an outstanding sample of intermediate type 
of stellar systems, specifically compact young star clus- 
ters, which are not available in the Galaxy. In such clus- 
ters primordial mass segregation can be suffi ciently observed 
from ground-based observations (e.g., Kontizas et al.lll998|; 
[Santiago et ai]|200TI: iKerber et al.H 2002: Kumai-etal. 200^ 
and they are also suitable for more detailed studies of this 
phenomenon from space with the high -resolving efficiency 
of the Hubble Space Telescope (e.g. Fischer et al. 1998|; 
de Griis et al. 2002; Sirianni et al. 2002; Gouliermis et alT] 
2004; Kerb er & Sant iago 2006). Such studies contribute sig- 
nificantly to the debate about the phenomenon of primordial 
mass segr egation and its orig in which which is still ongoing 
(see, e.g., 'Kissler-Patig et al.' 2007). It is interesting to note 
that while the existence of the phenomenon itself is being 
questioned (Ascenso et al. 2009), new methods are still being 
designed for its compre hensive identification and quantifica- 
tion (e.g.. lAmson et al.l l2b09a). 

In the first part of our study of stellar stratification we de- 
veloped and tested a new robust diagnostic method aimed at 
establishing the presence (or not) of this phenomenon in a 
cluster (Gouliermis et al. 2009, from here on Paper I). This 
method is based on the calculation of the mean-square radius, 
the so-called effective radius, of the stars in a cluster for differ- 
ent magnitude ranges, and the investigation of the dependence 
of the effective radii of each stellar species on magnitude as 
indication of stellar stratification in the cluster. 

1.1. The Effective Radius Method 

Stellar stratification affects the presently observable physi- 
cal properties of the star cluster, and thus diagnostic tools for 
the identification of the phenomenon are based on the detec- 
tion of mass- or brightness-dependent changes of such prop- 
erties. Previously established methods involve (i) the study 
of the projected stellar density distribution of stars of dif- 
ferent magnitudes; here, the differences in the exponents of 
the slopes fitted are seen as proof of mass segregation (e.g., 
ISubramaniam et all 119931) , (ii) the radial dependence of the 
Luminosity or Mass Function (LF, MF) slope; a gradient of 
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Fig. 1. — Distribution of tlie uncertainties of our pliotometry for all stars 
detected in both the F555W and F814W bands in the areas of the clusters 
of our sample. Uncertainties in F555W are plotted with filled black and in 
F814W in open grey circular symbols. 

this slope outwar ds from the cent er of the cluster indicates 
stratification (e.g.. lKing et alJ ll995). and (iii) the variance of 
the core radius of the cluster as estimated for stars in spe- 
cific magnitude (or mass) groups; a trend of this variance with 
brigh tness (or mass) signifies stratification (e.g..' Brandl et alj 
119961) . While these methods can provide useful information 
about the segregation in a cluster, all of them are sensitive to 
model assumptions, since they strongly depend on data fitting 
for parametrizing stellar stratification; i.e., functional expres- 
sions for the surface density profiles, or the MFs and LFs are 
required for the derivation of their slopes and core radii. 

Our method, developed in Paper I for assessment of stellar 
stratification in star clusters, is based on the notion of the dy- 
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namically stable Spitzer radius of a star cluster, defined as the 
mean-sq uare distance o f the stars from the center of the clus- 
ter (e.g.. ISpitzeiill958 h. The observable counterpart of this 
radius, the effective radius, is given by the expression: 



N 



(1) 



where r, is the projected radial distance of the /th stellar mem- 
ber of the cluster in a specific brightness range and the 
corresponding total number of stars in the same brightness 
range. Different calculations of the effective radius are per- 
formed, each for stars in different brightness ranges that cover 
the whole observed stellar luminosity function of the cluster. 
The simple application of this method can yield direct infor- 
mation about the spatial extent (radial distribution) of stellar 
groups in different magnitude (mass) ranges. In a star clus- 
ter, where stratification occurs, the segregated brighter stars 
are expected to be more centrally concentrated and the cor- 
responding effective radius should be shorter than that of the 
non-segregated fainter stars. As a consequence, stellar strat- 
ification can be observed from the dependence of the effec- 
tive radius of stars in specific magnitude (mass) ranges on the 
corresponding mean magnitude (mass). This method is suf- 
ficiently described and tested in Paper I, where we show that 
it performs efficiently in the detection of stellar stratification, 
provided that the incompleteness of the photometry has been 
accurately measured and the contamination by the field pop- 
ulation has been thoroughly removed. Our diagnosis method 
is also independent of any model or theoretical prediction, in 
contrast to those used thus far for the detection of mass segre- 
gation. 

In this second part of our study of stellar stratification we 
apply the effective radius method on deep imaging data ob- 
tained with the Wide-Field Channel (WFC) of the Advanced 
Camera for Surveys (ACS) onboard the Hubble Space Tele- 
scope (HST) of three young star clusters in the Large Mag- 
ellanic Cloud (LMC) for the investigation of primordial mass 
segregation in them. Young clusters in the Magellanic Clouds 
are quite different to each other, so that their study requires 
the detailed treatment of each cluster individually (see, e.g., 
§§[3l|4li, and therefore a statistical investigation of the phe- 
nomenon of mass segregation in a large sample of such clus- 
ters is not yet meaningful. The purpose of the present study is 
the detailed investigation with the use of a consistent method 
of individual selected clusters, in order to establish a compre- 
hensive comparative framework for the assessment of mass 
segregation in young LMC clusters. The objects of interest, 
clusters NGC 1983, NGC 2002 and NGC 2010, are selected 
among others also observed with ACS due to their youthful- 
ness and the variety in their characteristics. This selection is 
made in order to include in this comparative study clusters, 
which are very different from each other. In §|2]we describe 
the considered sample of LMC clusters and the corresponding 
data sets, as well as their reduction and photometry. We study 
the dynamical status of the clusters and discuss their structural 
parameters in §[3] and in §|4]we present the observed stellar 
populations, the decontamination of the stellar samples from 
the contribution of the local background field of the LMC and 
the stellar content of the clusters. The application of the diag- 
nostic method is performed for all three clusters in §|5] where 
we also develop a comparative scheme for the quantification 
of stellar stratification among different clusters, and we dis- 
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Fig. 2.— Completeness functions for NGC 1983 (top), NGC 2002 (mid- 
dle), and NGC 2010 (bottom) for both filters F555W (solid line) and F814W 
(dashed line). Average completeness functions over the whole observed fields 
of view are shown. It should be noted that completeness is a function of 
cluster-centric distance, as stars closer to the center of a cluster suffer from 
high crowding than those farther away. 

CUSS our results. We conclude on them in §|6l 

2. THE STAR CLUSTER SAMPLE 

The present study deals with stellar stratification in rich 
LMC star clusters. Such clusters are known to differ from 
each other in terms of, e.g ., their structural paramete rs, ages, 
luminosities and masses dMackev & Gilmorell2003h . since 
we focus on the phenomenon of primordial mass segrega- 
tion we limit our sample to young LMC clusters, and in or- 
der to use the most complete available data we select clus- 
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TABLE 1 

ACS/WFC OBSERVATIONS OF THE CLUSTERS OF OUR SAMPLE (HST PROGRAM 9891). 



Cluster 


R.A. 


DEC 


Filter 


Dataset 


Exposure 


Date 




(J2000.0) 


(J2000.0) 




filename 


time (s) 




NGC 1983 


05'"27""44?92 


-68°59'07"0 


F555W 


j8ne68req 


20 


October 7, 2003 








F814W 


j8ne68riq 


20 


October 7, 2003 


NGC 2002 


05'"30"'20?80 


-66°55'02"3 


F555W 


j8ne70b5q 


20 


August 23, 2003 








F814W 


j8ne70b7s 


20 


August 23, 2003 


NGC 2010 


05'"30"34?89 


-70° 49' 08 "3 


F555W 


j8ne71r5q 


20 


October 7, 2003 








F814W 


j8ne71r9q 


20 


October 7, 2003 



ters observed with ACS/WFC. The advantages introduced to 
crowded-field photometry of MCs star clusters by the use 
of ACS, and specifically its Wide-Field Channel, have been 
documente d extensively by studies based on such observa- 
tions (e.g., iMackev et al]|2006t iRochau et all2007l : iXin et alj 
120081) . Consequently, we select three young LMC clus- 
ters, NGC 1983, NGC 2002 and NGC 2010, observed with 
ACSAVFC. These clusters are different from each other, as far 
as their appearance, structure and stellar content is concerned. 
Considering that we also wish to investigate the identification 
and quantification of the phenomenon of stratification in var- 
ious cluster-hosting environments, we based the selection of 
these specific clusters on the on the differences between their 
local surrounding field star populations. We discuss in de- 
tail the differences of the clusters themselves in their structure 
and dynamical behavior in § |3] and in their stellar content in 
§ m The decontamination of the true cluster populations from 
the contribution of the general background field of the LMC, 
which determines the actual stellar content of the clusters and 
the stellar sample to be used in our method for assessment of 
stellar stratification in the clusters is performed in § 14.31 

2.1. Data Reduction and Photometry 

Our observations come from HST program 9891, a snap- 
shot survey of ^ 50 rich star clusters in the Magellanic 
Clouds. Imaging was obtained in Cycle 12 using ACS/WFC. 
As snapshot targets, each cluster was observed for one orbit 
only, resulting in a single exposure through the F555W filter 
and another through the F814W filter. Details of the observa- 
tions for the clusters considered in this paper may be found in 
Table □ 

The ACS WFC consists of a mosaic of two 4096 x 
2048 pixel CCDs covering an area of approximately 202" x 
202" at a plate scale of 0.05 arcsec per pixel. The two CCDs 
are separated by a gap roughly 50 pixels wide. The core of 
each cluster was placed at the center of WFC chip 1. This 
ensured that the inter-chip gap did not interfere with the main 
part of the cluster, and meant that each image reached a max- 
imum radius of approximately ^ 150" from the cluster cen- 
ter A small dither was made between the two exposures of a 
given target, to help facilitate the removal of cosmic rays and 
hot pixels; however with only two images per cluster it was 
not possible to eliminate the gap in spatial coverage due to 
the inter-chip separation. 

The data products produced by the standard STScI reduc- 
tion pipeline, which we retrieved from the public archive, 
have had bias and dark-current frames subtracted and are di- 
vided by a flat-field image. In addition, known hot-pixels and 
other defects are masked, and the photometric keywords in 
the image headers are calculated. We also obtained distortion- 
corrected (drizzled) images from the archive, produced using 



the PYRAF task MULTIDRIZZLE. 

We used the DOLPHOT photometry software (e.g.. lDolphinl 
l2000h . specifically the ACS modul^l to photometer the 
flatfielded (but not drizzled) F555W and F814W images. 
DOLPHOT performs point-spread function (PSF) fitting using 
PSFs especially tailored to the ACS camera. Before perform- 
ing the photometry, we first prepared the images using the 
DOLPHOT packages ACSMASK and SPLITGROUPS. Respec- 
tively, these two packages apply the image defect mask and 
then split the multi-image STScI FITS files into a single FITS 
file per chip. We then used the main DOLPHOT routine to 
simultaneously make photometric measurements on the pre- 
processed images, relative to the coordinate system of the 
drizzled F814W image. We chose to fit the sky locally around 
each detected source (important due to the crowded nature of 
the targets), and keep only objects with a signal greater than 
10 times the standard deviation of the background. The out- 
put photometry from DOLPHOT is on the calibrated VEGA- 
MAG scale of Sirianni et al. (2005), and corrected for charge- 
transfer efficiency (CTE) degradation. 

To obtain a clean list of stellar detections with high quality 
photometry, we applied a filter employing the sharpness and 
"crowding" parameters calculated by DOLPHOT. The sharp- 
ness is a measure of the broadness of a detected object rela- 
tive to the PSF - for a perfectly-fit star this parameter is zero, 
while it is negative for an object which is too sharp (perhaps 
a cosmic-ray) and positive for an object which is too broad 
(e.g., a background galaxy). The crowding parameter mea- 
sures how much brighter a detected object would have been 
measured had nearby objects not been fit simultaneously. We 
selected only objects with -0.15 < sharpness < 0.15 in both 
frames, and crowding < 0.25 mag in both frames. We also 
only kept objects classified by DOLPHOT as good star^ with 
formal errors from the PSF fitting less than 0.2 mag in both 
F555W and F814W. In Figure [T] typical uncertainties in the 
derived magnitudes of the stars detected in all three clusters 
are shown for both F555W and F814W filters. 

In order to assign a detection completeness fraction to each 
star, we used DOLPHOT to perform a large number of arti- 
ficial star tests. For each star in our final cleaned catalogue 
we first generated a list of 150 associated fake stars using the 
DOLPHOT utility ACSFAKELIST. We constrained these fake 
stars to lie close to the real star on the (F555W, F814W) color- 
magnitude diagram (CMD) - the fake stars were randomly 
and uniformly distributed in a box of size 2Av x 2Ac cen- 
tered on the real star, where Ay was the maximum of 0. 1 mag 

' The ACS mod ule of DOLPHO T is an adaptation of the photometry 
package HST PHOT 4Dolpl^ ^j2QQfl ). The package can be retrieved from 
http : / /pur cell . as ■ arizona . edu/ dolphot / 

^ Which have object type 1, as opposed to elongated or extended objects 
which have object types > 1 . 
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Fig. 3.— The ob.served ACSAVFC field.s-of-view of the clusters (a) NGC 1983, (b) NGC 2002 and (c) NGC 2010 (left panel), with the corresponding stellar 
density maps constructed with star counts on the photometric catalogs (right panel). Symbol sizes correspond to the brightness of the marked stars. These maps 
demonstrate that the selected clusters show a variety in size, structure, stellar density and ambient local field (see text in §[3). 



and 3 times the uncertainty in the nisss magnitude of the real 
star, and Ac was the maximum of 0.1 mag and 3 times the 
uncertainty in the (m^ss ~ iii&h) color of the real star 

We were also careful to constrain the spatial positions of the 
fake stars. Often this is done so that the fake stars lie at simi- 
lar cluster-centric radii to the real star under consideration, but 
with unconstrained position angle. This saves on computation 
time because the calculations are then valid for groups of real 
stars at a time, rather than simply on an individual basis. How- 
ever, because we are dealing with very young clusters in the 



present work, this implicit assumption of spherical symmetry 
is not necessarily valid. These objects could potentially still 
be embedded in gas or dust and thus exhibit a strongly vary- 
ing background. They are also not very dynamically evolved, 
so may yet possess azimuthal asymmetries in stellar density 
as well as numerous randomly placed bright stars. Therefore, 
we adopted the more computationally expensive method of 
generating the list of fake stars locally around each real star 
on the CCD - in this case randomly and uniformly within a 
radius of 5". In doing this we were careful to account for the 
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TABLE 2 

Characteristics of the clusters in our 

SAMPLE derived FROM OUR PHOTOMETRY. 



Cluster 


Size 




Age 


Av 




(pc) 


(stars arcsec ^) 


(Myr) 


(mag) 


NGC 1983 


7.3 


0.72 


28 


0.17 


NGC 2002 


9.4 


0.65 


18 


0.19 


NGC 2010 


12.5 


0.54 


159 


0.17 



Note. — Sizes are derived from circular annuli en- 
compassing the 3(7 isopleth in the iso-density contour 
maps of the clusters in the observed regions. Stellar den- 
sity /3o- is measured within these annuli and are given 
for comparison with the field stellar density of Table |3] 
Age and interstellar extinction for the clusters are esti- 
mated from evolutionary model fitting on the CMDs of 
stars within the of the clusters. Sizes are given in phys- 
ical units after a conversion 1' ~ 15 p c assuming a dis- 
tance modulus of (m-M)o :^ 18.5 mag iClementini et alj 
120031: IAlvesl200llSchaefeil2008ft for all clusters. 

edges of the chips, including the inter-chip gap. 

We then ran DOLPHOT in artificial stars mode on each list 
of fake stars. When operating in this way DOLPHOT takes the 
next fake star from the list, adds it to the input images, and 
then solves for the position and photometric properties of the 
star using exactly the same parameters as for the original (real 
star) analysis. Once complete, we took the output photometry 
and ran it through our quality filters, just as we did for the list 
of real stars. Finally, we calculated the completeness fraction 
for a given real star by determining how many of its associated 
fake stars had been successfully 'found'. We considered a 
fake star to be found if it was recovered in both F555W and 
F8 14W and passed successfully through our quality filter, and 
provided its output position lay within 2 pixels of its input 
position. Figure |2] presents the completeness function of our 
photometry for all considered clusters in both F555W (solid 
Une) and F814W (dashed line) filters. 

3. CLUSTER MORPHOLOGIES AND STRUCTURES 

3.1. Morphology of the Clusters 

The stellar charts of the observed ACSAVFC fields of the 
clusters are shown in Figure |3] (left panel). Our observations 
show clear differences in the appearance of the three clusters. 
In order to reveal the spatial distribution of the detected stars 
and the morphology of the clusters, we performed star counts 
on our photometric catalogs. The constructed 'iso-density' 
contour maps for each region are also shown in Figure|3](right 
panel). The coordinate system used for these maps is that 
of the original ACSAVFC pixel coordinates. The star counts 
were performed in a quadrilateral grid divided in elements 
with sizes ^ 190 x 190 WFC pix^ each, which correspond to 
about 10" X 10" (or « 2.5 pc x 2.5 pc at the distance of the 
LMC). The first isopleth in the contour map indicates the level 
of Icr above the mean background density, with a being the 
standard deviation of the background density. The subsequent 
levels are drawn in steps of la. All contours with density 
equal or higher than the 3<t level, considered as the statisti- 
cally significant density threshold, are drawn with thick lines. 
All three star clusters are revealed in these maps as the dom- 
inant stellar concentrations in the observed regions. These 
maps demonstrate that while NGC 1983 is a centrally concen- 
trated compact cluster, NGC 2010 is a quite large and loose 
cluster, with NGC 2002 being an intermediate case of a large 
compact concentration. Moreover, NGC 2002 appears to be 
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Fig. 4.— Radial suiiace density profiles of NGC 1983 (top) NGC 2002 
(middle), and NGC 2010 (bottom) with the best-fitting King and EFF mod- 
els superimposed. The en'or bars represent Poisson statistics. Left. Density 
profiles of all detected stars, including the local field of the LMC. The best- 
fitting EFF model is over-plotted with a solid blue line. Right. Stellar surface 
density profiles after subtracting the background density level estimated from 
the best-fitting EFF model shown in the plots on the left. Both EFF and King 
models are overlaid. Since these clusters are not tidally truncated, the King 
models are used only tentatively for an estimation of a limiting radius for the 
clusters. Both models, however, seem to fit very well the density profiles at 
the inner parts of the clusters deriving comparable core radii. The derived 
structural parameters based on the best-fitting EFF (solid blue line) and King 
(dashed green line) models are given in Tabled 

the most spherical cluster in the sample, since NGC 1983 is 
more elUptical and NGC 2010 is more amorphous, for which. 
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as indicated by the isopleths of density > 3fT, its most com- 
pact part is off-center. The sizes of the clusters as derived 
from circular annuli considered to encircle the 3(t isopleth for 
each cluster are given in Table |2] The corresponding surface 
stellar density of each cluster is also given in this Table. 

3.2. Stellar Surface Density Profiles 

The observed fields cover quite large areas around the clus- 
ters in our sample, allowing the construction of the stellar sur- 
face density profiles of the clusters at relatively large distances 
from their centers. We divided the area of each cluster in con- 
centric annuli of increasing steps outwards and we counted 
the number of stars within each annulus. We corrected the 
counted stellar numbers for incompleteness according to the 
corresponding completeness factors on a star-by-star basis. 
Completeness is a function of both distance from the center of 
the cluster and magnitude. Let A^, ^ be the completeness cor- 
rected stellar number within the ith annulus. We obtained the 
stellar surface density, /), by normalizing this number to the 
area of the corresponding annulus as /, = A^,.c/A,, with A, be- 
ing the area of the annulus. The observed fields do not always 
fully cover the complete extent of the clusters, and, therefore, 
for the truncated annuli, we considered only the available area 
for the estimation of the corresponding surface density. 

The "raw" radial stellar density profiles, meaning the sur- 
face stellar density as a function of distance from the cen- 
ter of each cluster, f{r), show a smooth drop outwards away 
from the centers of the clusters. They drop to a uniform level, 
which represents the stellar density of the LMC field in the 
vicinity of the clusters, and it is measured by fitting the models 
of Elson, Fall & Freeman ( 1987) to the stellar surface density 
profiles, log /(log r), (FigurelH left) as described below. In the 
case of NGC 2002 there is a small increase in the outskirts of 
the profile due to another neighboring small stellar concentra- 
tion (see Figure [3]). We do not consider this increase in our 
subsequent analysis on the structure of this cluster The errors 
in the profiles reflect the counting uncertainties, representing 
the Poisson statistics. 

3.3. Structural Parameters 



We applied both the empirical mod e l by [Kind (119621) and 
the model of lElson. Fall & Freeman I (Il987h (from here-on 
EFF) to the stellar surface density profiles of the clusters in 
order to obtain their structural parameters. The EFF model is 
suited for clusters which are not tidally limited, while King's 
empirical model represents tidally truncated clusters. Both 
models provide the opportunity to derive accurate character- 
istic radii for the clusters. The core radius (r^) describes the 
distance from the center of the cluster where the stellar density 
drops to the half of its central value and the tidal radius (rt) is 
the limit where the stellar density of the cluster drops to zero. 
The best-fitting EFF and King profiles are found by perform- 
ing the Levenberg-Marquardt least-squares fit to the consid- 
ered functions. This fit was performed with IDLQ (Interactive 
Data Language) with the use of the specially developed pro- 
cedure MPFITFUN (Markwardt 2008). For the application of 
King's model the stellar surface density of the cluster alone 
(with no contamination from the field) is necessary, while the 
EFF model does not require any field subtraction. We first 
apply the latter in order to estimate the core radius and the 
uniform background level for each cluster The subtraction of 

^ |http : / /www ■ it t vis . com/Product Service s/IDL . aspx] 



this density level from the measured surface density at each 
annulus gives the surface density profile of the cluster alone, 
from which we derive the core and tidal radii. 

3.3.1. Best Fitting EFF Profiles 

From studies of young LMC clusters by EFF it appears that 
these clusters are not tidally truncated. These authors devel- 
oped a model more suitable to describe the stellar surface den- 
sity profile of such clusters: 



/(r) = /o(l + rVa2p/2^ /field, 



(2) 



where /o is the central stellar surface density, a is a measure of 
the core radius and 7 is the power-law slope which describes 
the decrease of surface density of the cluster at large radii; 
f(r) cx r"'''/^ for r^a. The uniform background density level 
is represented in Eq. (|2| by /field. The measured parameters a, 
7 and /field, derived from the fitting procedure on our data are 
given for each cluster in Table[3] According to EFF model, the 
core radius, r^, is given from Eq. (|2]i assuming no contribution 
from the field as: 



rc = ail 



2/7 _ 1)1/2 



(3) 



The estimated core radii of the clusters derived from Eq. (|3]l 
are also given in Table |3] The best-fitting EFF models are 
shown superimposed (blue lines) on the stellar surface density 
profiles of Figured 

3.3.2. Best fitting King Profiles 

In order to construct the density profiles of the clusters 
alone with no contribution from the field we subtract the back- 
ground density level, /geid, from the raw profiles of Figure]?] 
(left). We then use the field-subtracted profile to deri ve an 
indicative tidal radius of the clusters as described by I Kind 
( 119621) . According to this model the density profile of a tidally 
truncated cluster is given as: 



/(r)oc f 

Hl + (^ 



[l + (^)2]2 [l+(i)2]2 



(4) 



where / is the stellar surface density, r^ and r^ the core and 
tidal radius, respectively and r is the distance from the center. 
The tidal radius, which represents the outskirts of the cluster, 
is found from the formula 



m = Ml/r-l/r,f 



(5) 



and the core radius, which describes the inner region of the 
cluster, is given from 



/('•) = 



./o 



l+(r/re)2' 



(6) 



/o describes again the central surface density of the cluster and 
/i is a constant. The best-fitting King profiles apart from core 
and tidal radii for the clusters, deliver also the concentration 
parameters, c, defined as the logarithmic ratio of tidal to core 
radius c = log(rt/rc), which refers to the compactness of the 
cluster The derived rt and c for our clusters are given in Ta- 
ble [3] while derived from the best-fitting King profiles are 
omitted, as they are found in excellent agreement with those 
derived above from the EFF profiles. The best-fitting King 
models of our clusters are shown with green dashed lines su- 
perimposed on the radial surface density profiles of Figure ]4] 
(right). 
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TABLE 3 

Results of the application of both EPF and King models fit to the observed radial surface density profiles of the clusters in 

OUR sample. 









EFF Profiles 






King Profiles 


Cluster 


a 


7 


/field 




n 


c 




(pc) 




(stars arcsec"^) 


(pc) 


(pc) 


logrt/rc 


NGC 1983 


1.68 ±0.28 


1.69 ±0.32 


0.23 ± 0.02 


1.89 ±0.43 


21.69 


0.8 


NGC 2002 


3.92 ± 0.68 


3.91 ± 1.00 


0.09 ± 0.01 


2.56 ± 0.61 


22.33 


0.9 


NGC 2010 


13.32 ±6.03 


9.81 ± 8.04 


0.21 ±0.01 


5.19 ±3.34 


28.13 


0.7 



Note. — Parameters a and 7 are derived from the best-fitting EFF profile according to Eq. J2j. Core radii, r^, are derived from the same models with 
Eq. (5). Tidal radii, rt, are measured with Eq. l|5lfrom the best-fitting King models. Concentration parameters, c, are estimated from the measured rt 
and core radii derived from King models and Eq. l|6). All radii are given in physical units after a conversion 1' ~ 15 pc assuming a distance modulus of 
(m-M)o ~ 18.5 mag iClementini et al. 2003; Alves 2004; Schaefer 2008) for all clusters. 



In general, both King and EFF models are in good agree- 
ment to each other in the inner and intermediate regions of 
the clusters, which are equally well fitted by both models. 
However, based on the fact that young compact clusters in the 
LMC are not tidally truncated, the EFF model should be con- 
sidered as the best representative for them. On the other hand, 
the decontamination of the observed stellar samples from the 
populations of the local LMC field requires t he d etermination 
of a limiting radius, Rw^, for the clusters (§ I4.3l l. Since EFF 
models do not predict such a radius, we treat the tidal radii 
derived from the application of King's theory as this radius 
for each cluster. It should be noted, though, that our clusters 
are not expected to be dynamically relaxed and therefore, the 
results on their tidal radii based on the application of King 
profile fitting should not be taken literally but only tentatively 
for the determination of an indicative Jimit radius for the clus- 
ters, and for reasons of comparison. For example, as seen in 
the charts of Figure [3] NGC 1983 cannot be directly distin- 
guished from the rich stellar field in which it is embedded, as 
opposed, e.g., to NGC 2002, which shows up as a compact 
concentration easily separated from its surroundings. 

4. STELLAR CONTENT OF THE CLUSTERS 

In the stellar charts of the three observed regions, shown 
in Figure [3] (left), the star clusters immediately appear as the 
most prominent stellar concentrations. These maps are con- 
structed only from stars detected with photometric uncertain- 
ties (T < 0.1 mag in both F555W and F814W filters. In our 
subsequent analysis we consider only these stars for each ob- 
sei-ved field. The ACSAVFC field-of-view of NGC 1983 cov- 
ers 9 641 stars with good photometry, the one of NGC 2002 
4 340, and that of NGC 2010 9 006 such stars. 

4.1. Color-Magnitude Diagrams 

The m555 -OT8i4, m^w Color-Magnitude Diagrams (CMDs) 
of the total observed samples of stars for all three clusters are 
shown in Figure |5] (top panel). These CMDs of the VVTioJe 
Areas demonstrate the existence of a significant mixture of 
different stellar populations in each region; from young stars, 
which cover the main sequence (MS) to the turn-off, belong- 
ing to the clusters and their surrounding field, to the evolved 
Red Giant Branch (RGB) and Red Clump (RC) stars, repre- 
senting mostly the local background LMC field. The CMDs 
of Figure |5] indicate that our photometry provides accurate 
magnitudes for stars with msi4 ^ 22 mag. It is interesting 
to note that as seen from the CMDs of Figure |5] (top), as well 
as from the stellar charts of Figure [3] the general field of the 
LMC in the observed regions shows to vary in stellar density. 
Specifically, both NGC 1983 and NGC 2010 are located in 



much denser stellar fields than NGC 2002, with the local field 
of NGC 1983 being particularly rich in red field stars. This 
is akeady demonstrated by the measurements of the surface 
s tellar density, /field, of the general field of all three clusters in 
§[33](Table[^ column 4). 

The tidal radii, rt, of the clusters, as they are derived in 
§ 13.31 basically define the spatial extent of each cluster, and 
therefore we treat the stars that are located at radial distances 
from the center of the clusters larger than r^ as the best rep- 
resentative stellar populations of the surrounding fields of the 
clusters. The corresponding CMDs, for which we use the term 
Field CMDs, are shown in Figure |5] (middle panel). In these 
CMDs the features of both RGB and RC of the general field 
of the LMC are well apparent in all three observed regions. 
These populations are well documented as members of the 
LMC field in various HST s tudies of the star forma tion history 
(SFH) of this ga laxy (e.g., 'Geha et al.in" 998.: .Holtzman et alj 
[l999l : ICastro"eta l. 2001; Smecker-Hane eLdJ|20Q2|). Specif- 
ically, the SFH of the LMC is found to be continuous with 
almost constant star formation rate (SFR) for over ^10 Gyr, 
with an i ncrease of the SFR at around 1-4 Gyr ago (for a sum- 
mary see i Gouliermis et al. 2006. their § 3.1). This old LMC 
field population is present in all the Field CMDs of Figure |5] 
Since it covers the same parts in all three of them, we define 
the age-limits of these stars only on the CMD of the field of 
NGC 2010. With the use of the Padova grid of evolutionary 
models (Girardi et al. 2002) for insignificant extinction and 
metallicity and distance typical for the LMC these limits are 
found to be 1 ^r/Gyr^ 10. 

Moreover, the Field CMDs of Figure |5] include also young 
MS stars, clearly suggesting that the general regions, where 
all three clusters are located, include young stellar popula- 
tions that are not members of the clusters. This suggestion 
does not contradict previous studies on such clusters and their 
surroundings. Indeed, young compact clusters in the LMC 
are known to form together in larger structu res of young stel - 
lar populations, the Stellar Aggregates (e.g.. lOev et al.ll2008h . 
the size of which seems to cor relate with the duration of star 
formation in them (e.g., E fremov & Elmegreenlll999l) . Such 
clusters are not isolated a s they may be form ed in binary or 
multiple cluster systems ( Dieball et alJl2002h . An inspection 
of the general regions where the clusters of our sample are lo- 
cated verifies that none of them is an exception to this general 
rule. Wide-field images of the regions around the clusters and 
the whole extent of the LMC in B, R and / were retriev ed from 
the SuperCOSMOS Sky Survey ("Hamblv et al .Il2001h . avail- 
able at the Aladin Sky Atlas (Bonnarel et a rT l2000h . These 
images show that indeed all three clusters are members of 
larger young structures. Specifically, NGC 1983 is located 
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Fig. 5. — The 111555- ™iu vs. m8i4 CMDs of the stars detected with the best photometry (uncertainties cr < 0.1 mag in both filters) in the observed HST 
ACSAVFC areas of the LMC clusters of our sample: NGC 1983 (left), NGC 2002 (middle) and NGC 2010 (right). Top: The complete CMDs of all stars detected 
in the whole observed areas, which demonstrate the mixture of the young stellar populations of the clusters (MS) with those of their surrounding field, and 
the evolved ones of the general LMC field (RGB, RC). Middle: The CMDs of the stars located outside the of the clusters, representing the most probable 
populations of the local ambient of each cluster. Bottom: The CMDs of th e stars com prised within the of each cluster, which represent best the true CMDs of 
t he cl usters. Isochrones from the grid of evolutionary models bv .Girardi et alj 120021) ai'e overlaid for an estimate of the age of the clusters and the old field (see 

§E3)- 



at the c entral par t of the star forming Shapley Constellation 
JJ (IShaplevJ95illMcKibben Nail & Shaplev 1953) that coin- 
cides with the supergiant shell (SGS) LMC 3 (Meaburn 19801) . 
NGC 2002 belongs to the northern part of Constellation 777 or 
SGS LMC 4, and NGC 2010 is located in one of the eastern 
nebular filaments of SGS LMC 9 in Constellation IX. 

In order to compare these MS field populations with the 
actual stellar populations of the clusters we select the stars 
contained in the core radii of the clusters (see § 13.31 Table [3] 
column 5), as the most probable stellar members, and we 
construct the corresponding CMDs, which are shown in Fig- 
ure |5] (bottom panel). A comparison of the Field CMDs with 
the CMDs of the clusters within their core radii shows that 
there is a confusion between the young field populations and 
those of the clusters, especially for NGC 2002 and NGC 1983. 
The case of NGC 2010 is more straightforward, as it can be 
seen that while the cluster is somewhat evolved (see discus- 



sion below), it seems that it is 'embedded' in a younger field 
with brighter MS stai's. On the other hand, NGC 1983 and 
NGC 2002 are more complicated cases as both field and clus- 
ters seem to comprise the same type of MS stars. Taking into 
account the similarities of field stars to cluster members in the 
CMDs, the use of our method for the statistical decontamina- 
tion of the observed stellar samples from the field populations, 
discussed in § 14.31 becomes quite important. As we show in 
§ l4.4l this method allows us an accurate statistical determina- 
tion of the most probable stellar members of each cluster. 

4.2. Ages of the Clusters 

We use the CMDs within the core radii of the clusters to de- 
rive an estimation of their ages. In the corresponding CMDs 
of Figure |5] (bottom panel), the best-fitti ng isochrones from 
the Padova grid of evolutionary models dGirardi et al.ll2002h 
are also plotted. The process of isochrone fitting of the ob- 
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served CMD results in an accurate measurement of the age 
of the system, providing that there are distinguished features 
in the CMD to be fitted. Moreover, a good estimate of the 
distance and the interstellar reddening of the cluster is re- 
quired. For clusters in the LMC, visual extinction is known 
to vary around very low values (e.g., de Griis et al. 2002; 
iGouliermis et al. 1 120041 : iKerber & Santia go 2006), wh ile the 
distance of the gala xy is also well defined (e.g., lAlvesI 
l2004l : lSchaefeai2{)08h . Indeed, for our age estimation we 
had to vary little both of these parameters around typi- 
cal values. Since the ACS photometric system is slightly 
different from the standard Johnson-Cousins system, we 
estimate the mean extinction in the F555W band, A555, 
which we arbitrarily symbolize as Ay, u s ing th e compu- 
tations by iDa Rio. Gouliermis & Henn ing* ('2009'). These 
authors derive a conversion R^ss = A^ss/Eim^ss - I'nau) — 
2.18 and A5s- i /Asi4 1.85 assuming the extinction law of 
iCardelh et alf (11981 with Ry ~ 3.1. In all cases we assume 
the typical metall icity of the LMC of Z ~ 0.3 - 0.5 Zq (e.g., 
IWesterlundll 1 997h . and therefore we use the models designed 
for Z = 0.008. The derived ages and indicative visual extinc- 
tion for all three clusters are given in Table |2l Let us consider 
the results for each cluster individually: 

(i) NGC 1983: The CMD of the most prominent member 
stars of the cluster consists only of the MS, indicating that 
the cluster is quite young, but making the accurate determi- 
nation of its age difficult. Indicatively, we apply the youngest 
isochrone that fits the brighter MS stars with ;«8i4 ^ 14 mag, 
which shows that the cluster has an age about 28 Myr The 
selected isochrone fits best the blue part of the MS for a 
nominal distance modulus (m-M)o ~ 18.45 mag and red- 
dening of ^(msss -OT814) ^ 0.05 mag, which corresponds to 
Ay — 0.17 mag. 

(ii) NGC 2002: The CMD of NGC 2002 within the of 
the cluster includes the significant number of five Red Su- 
pergiants (RSG) located at 2.2 ^ (^555 - mgi4) ^ 2.6 mag 
and 11 ^ m8i4 ^ 12 mag, which helps to constrain ac- 
curately the age of the system. However, the isochrone for 
metallicity Z = 0.008 that seems to fit best this cluster with 
an age of r ~ 18 Myr, cannot reach the RSG stars, which 
thus require the consideration of models of higher metallic- 
ity. Indeed the next available models from the Padova grid 
are those of almost solar metallicity of Z = 0.019. The cor- 
responding isochrone for the same age (plotted with a green 
line in Figure |5]l does fit the RSG, implying that either there 
is a significant change in the metallicity of the cluster, or that 
the evolutionary models need to be fine-tuned for RSGs in 
low metallicities. In any case, our data give a fixed age for 
NGC 2002 of ai-ound 18 Myr. It should be noted that while 
the high metallicity model requires a modest Ay ~ 0.12 mag 
to fit the CMD, the low metallicity model predicts a higher 
extinction of Ay — 0.19 mag. 

(iii) NGC 2010: The CMD of Figure |5] (bottom) for 
NGC 2010 clearly suggests that this is the oldest cluster 
in our sample. The existence within of bright stars of 
m8i4 ~ 15 mag, evolved away from the MS allows us an ac- 
curate determination of the cluster age of r ~ 159 Myr for 
(m-M\) ~ 18.5 mag and Ay ~ 0.17 mag. However, the ap- 
pearance in this spatially constrained CMD of a prominent 
RGB population, clearly suggests that the cluster may com- 
prise a mixture of different stellar populations. 

4.3. Field Decontamination 



As we discuss earlier and seen from the CMDs of Fig- 
ure |5] all three observed regions suffer from significant con- 
tamination by field stars. With no information about cluster- 
membership probabilities for the observed stars, obtained, 
e.g., from radial velocities and/or proper motions, the quan- 
titative decontamination of the clusters from the field stars on 
a statistical basis becomes a fundamental method to obtain the 
CMD of the complete sample of true stellar members of each 
cluster We decontaminate, thus, the observed CMDs from 
the contribution of the stellar members of the surrounding lo- 
cal field of the LMC with the use of a sophisticated random- 
subtraction technique. 

This field-subtractio n algorithm, originally designed by 
iBonatto & Bical (120071) . is already developed in Paper I, where 
it is also thoroughly described. According to this technique, 
the successful decontamination of the clusters from the field 
is based on the definition of the limiting radius, Rxim, for each 
cluster For the clusters in our sa mple, /?iim is given by the 
tidal radii, rt, estimated in § 13. 3. 21 with the use of King's the- 
ory. All stars located at distances r > from the center 
of the cluster are treated as field stars, while the rest (with 
f < ^lim) are considered as the most probable cluster-member 
stars. Assuming a homogeneous field-star distribution, the 
number density of field stars is calculated in the Field CMD 
(Figure|5l middle) constructed for stars with r > Rn^n, consid- 
ering also the observational uncertainties in the photometry. 

The decontamination process is then applied in three steps 
for each cluster: (i) Division of the CMDs of the stars with 
r < R]i^ (cluster region) and of those with r > R^^ (offset re- 
gion), respectively, in 2D cells with the same axes along the 
m8i4 and (mg55-m8i4) directions, (ii) Calculation of the ex- 
pected number density of field stars in each cell in the CMD 
of the offset region, (iii) Random subtraction of the expected 
number of field stars from each cell from the CMD of the 
cluster region. The algorithm is applied several times using 
different cell sizes (Am8i4, A(m555-m8i4)) so that many dif- 
ferent decontamination results are derived. From these results 
we calculate the final probability that a star is identified as a 
'true' cluster member, and the final decontaminated CMD of 
the cluster contains only the stars with the highest probability 
of being cluster members. This process allows us to minimize 
any artificial effects intrinsic to the method itself. A detailed 
description of the mathematical formulation of our method for 
field-subtraction is given in Paper I. 

4.4. Stellar Populations in the Clusters 

We performed the field-subtraction process, as described 
above, 15 times for each cluster considering combinations for 
the CMD cell sizes of Amgn = 0.5, 1.0, 1.5, 2.5 and 3.0 mag 
and A(m555-m8i4) = 1.0, 1.5 and 2.0 mag. We obtained, thus, 
the probability for each star in each sample of being a mem- 
ber of cluster population. Since our method for diagnosis 
of stellar stratification is completeness-sensitive in the sub- 
sequent analysis we limit the sample of stellar members of 
all three clusters to those detected with our photometry with 
completeness > 70%. The derived stellar samples comprise 
424 stars in NGC 1983, 1 102 stars in NGC 2002 and 2 265 
stars in NGC 2010 with probability larger than 70% of be- 
ing true cluster members. Figure |6] shows the correspond- 
ing field-star-decontaminated CMDs of the cluster regions for 
all clusters. Stars found by the field-subtraction method with 
probability 100% of being cluster members are plotted with 
thick points. From these CMDs it can be seen that our statis- 
tical subtraction of the CMD contamination by the field stars 
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Fig. 6. — The 1x1555- msw vs. m8i4 CMDs of the stars retained as possible members of the clusters after the statistical removal of the contribution of the LMC 
field from the originally observed CMDs of Figure|5](top panel). The Field CMDs of the offset regions around the clusters ( Figu re [3] middle panel) are used as 
representative of the contaminating populations. Stars with probabilities > 70% of being true members of the clusters (see § 14. 3t are considered for these plots. 
Stars identified with 100% membership are plotted with thick points. 



was efficient in removing both young field populations of the 
immediate surroundings of the clusters, and old populations 
of the general LMC field. Specifically, for NGC 1983 the MS 
stars remaining in the CMD of Figure |6l after the applica- 
tion of the field-subtraction, form a sharp narrow MS down to 
m8i4 ~ 21 mag. This MS resembles that of the CMD within 
the rc of the cluster shown in Figure |5] (bottom), with no indi- 
cation of the broad lower MS stars of the Field CMD of Fig- 
ure |5] (middle). The CMD of NGC 2002 also demonstrates 
a very sharp MS similar to the CMD of the cluster. The 
RSGs located within the of the cluster also remained in the 
CMD of Figure |6] with 100% of membership probability. As 
far as NGC 2010 is concerned, the statistical field subtraction 
r emov ed completely the bright MS young stars identified in 
§ 14.11 as young field stars, while the older MS and evolved 
stars remained with high probabilities of being true cluster 
members. 

It should be noted that the lower main sequence of the 
CMDs of the clusters may be also sensitive to any potential 
pre-main sequence (PMS) population that might exist in such 
young clusters. However, due to observational constrains, 
we would be able to detect only the turn-on of these PMS 
stars, which is known to coincide with the MSTO (see, e.g., 
iNota et al.1 [20061 : iGouhermis et al.l 12007 ). very close to our 
detection limit. As a consequence, any PMS population in 
the clusters would be seriously affected by the incomplete- 
ness in our photometry and therefore does not have any sig- 
nificance in the measured stellar numbers. Another bias that 
may affect the final numbers of low main sequence cluster- 
members is differential reddening that could displace stars to 
the red-faint part of the CMD. However, differential redden- 
ing in young LMC clusters i s known to account f or no more 
than A£(g-y) ~ 0.1 (e.g.. iDirsch etal . 2000; Baume et alj 
120071) . which is well covered by the color uncertainties in our 
photometry at the low main sequence. As a consequence, dif- 
ferential reddening is not considered in our treatment for field 
subtraction. It is interesting to note that while the Cluster 
CMDs of Figure for both NGC 2002 and NGC 2010 show 
the same features with the corresponding CMDs of the clus- 
ters (Figure |5] bottom), that of NGC 1983 shows clear indi- 
cations of multiple young populations with msu ^ 15 mag 



Indeed, as we discuss in § 14.11 multiplicity is a typical char- 
acteristic of young LMC clusters, and naturally it may lead to 
a merging process and mixin g of different stellar population s 
within one cluster (see e.g.. iPortegies Zwart & Ru^ 120071) . 



As we cannot verify or reject this scenario for NGC 1983, 
we treat all 424 stars that meet our membership criteria as 
true members of the cluster. Our diagnostic method for stellar 
stratification is applied in the following section on the final 
catalogs of stellar members of the clusters, with 70% mem- 
bership, the CMDs of which are shown in Figure |6l 

5. DIAGNOSIS OF STELLAR STRATIFICATION IN THE CLUSTERS 

The effective radius method for the diagnosis of stellar 
stratification is based on the assumption that the observable 
counterpart of the Spitzer radius, the effective radius, of stars 
in a specific magnitude (or mass) range will be a unique func- 
tion of this range if the system is segregated. In Paper I the 
comparison of the results of this method with those from the 
'classical' methods of the radial dependence of the clusters 
LFs and MFs, showed that the effective radius method be- 
haves more efficiently in the detection of the phenomenon, 
providing direct proof of stellar stratification in star clusters 
without considering any model or theoretical prediction, and 
with no application of any functional fit. In this section we ap- 
ply this method for the detection and quantification of stellar 
stratification to the young LMC clusters of our sample. 

5.1. Effective Radii of the Star Clusters 

We calculated the effective radii, r^a, of the stars in differ- 
ent magnitude ranges, according to Eq. ([U, using the com- 
pleteness corrected numbers of stars for all three clusters. We 
performed this calculation only for the star-members of the 
clusters, as they are established in § l4.4l (CMDs of Figure|6]), 
and we used their brightness measured in the F814W filter, 
since this filter provided the most complete and accurate pho- 
tometry (see § 12. 11 1. In order to identify any functional rela- 
tion between r^ff and the corresponding brightness range, we 
binned the stars according to their magnitudes. Since the mea- 
surement of Teff depends strongly to the number of stars per 
bin, this process naturally inherits systematic uncertainties to 
the derived relations, because of the fact that some stars could 
belong to neighboring bins due to photometric uncertainties. 
Therefore, we performed the same binning process per cluster 
several times by changing the bin sizes and/or the magnitude 
limits, within which the stars are binned. We found that a rea- 
sonable bin size, which provides good statistics and allows for 
the detailed observation of the dependence of ^eff on magni- 
tude, is that of 0.5 ^ Am8i4 ^ 1.0 mag. For this range of 
bin sizes we did not find any significant change in the general 
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Fig. 7. — Estimated effective radii of stars witliin different magnitude 
ranges vs. tlie corresponding mean magnitude for all three clusters. 

trend of reff(m8i4), which is also found to be independent of 
the selected magnitude limits for binning the stars. 

In Figure |7] we indicatively show the relations between 
the computed effective radii and the corresponding magni- 
tude bins derived from two different binning processes with 
bin sizes of 1 mag and magnitude ranges shifted by 0.5 mag 
from each other. The combination of both distributions with 
these specific binning characteristics (each plotted with differ- 
ent symbols) provides measurements of r^ff for every 0.5 mag 
range. The derived trends of reff(mgi4), shown in these plots, 
are not different from the corresponding distributions derived 
for Am8i4 = 0.5 mag, except of being somewhat smoother. In 
general the relations of Figure |7] are representative of those 
found from any applied binning, and therefore we consider 
them in the following section for the derivation of our results 



concerning stellar stratification in the clusters. 

5.2. Stellar Stratification in the Star Clusters 

The r^f{{m%n) graphs of Figure |7] are plotted on the same 
scale for reasons of comparison between the clusters. In these 
plots it is shown that the effective radius does behave as a 
function of magnitude, but not necessarily providing proof 
of stellar stratification for all clusters. A definite indication 
of stellar stratification would require the effective radii of 
the faintest stars to be systematically larger than those of the 
brightest ones for the whole observed brightness range. This 
is not the behavior we observe at least for two of the three 
clusters of our sample. Indeed, as we discuss in Paper 1, reft is 
not a monotonic function of brightness, and mass segregation, 
if any, is not expected to behave in the same manner in every 
cluster, but rather to depend on the intrinsic stellar character- 
istics of each of them. For example, while in NGC 2002 there 
is a trend of fainter stars having systematically larger effec- 
tive radii for the whole luminosity range, this is not the case 
for NGC 1983 that shows this trend for only a limited magni- 
tude range, and certainly not for NGC 2010, which one may 
argue is not even segregated. 

Specifically, for NGC 1983, the brightest magnitude bins, 
down to mgi4 ~ 16 mag, seem to include stars that are located 
away from the central part of the cluster. Indeed, the brightest 
stars in the Cluster CMD (Figure |6]l with m8i4 ^ 15.5 mag 
are found to be distributed within the largest r^fi of about 1 .'2 
away from the center of the cluster and not closer than about 
0f6. On the other hand, taking into account the remaining 
magnitude bins of > 15.5 mag, the graph of reff(m8i4) 
for NGC 1983 provides clear indications that the cluster is 
segregated with an almost monotonic dependency of larger 
Teff to fainter stars. 

The graph for NGC 2002 is quite different, giving evidence 
that the cluster is segregated in the whole extend of observed 
stellar luminosities. There is, however, a dependency of the 
appearance of a relation between reft and brightness on the 
considered magnitude range. Indeed, while the three brightest 
magnitude bins, comprising six stars (the five RSG and the 
brightest MS star), correspond to the smallest effective radii, 
the next fainter bins, which include six more MS stars, show 
that these stars are located well away from the center of the 
cluster. Moreover, from the plot of Figure |7] for the fainter 
stars, it is almost self-evident that they are segregated as the 
corresponding relation reff(mgi4) demonstrates an increase of 
Teff with brightness. 

The behavior of this relation for NGC 2010 is quite differ- 
ent from both the previous cases, as there is no general trend 
for larger effective radii with fainter magnitudes observable, 
and thus the cluster behaves as if there is no significant strat- 
ification occurring. Specifically, apart from the fact that the 
three brightest magnitude bins are inconclusive, for more faint 
stars there is a "plateau" in the function reff(m8i4), which ap- 
pears rather flat for magnitudes down to mgi4 2± 19 mag, while 
for even fainter stars, down to the faintest observed magni- 
tude, the dependence of reff over mg^ appears again, giving 
evidence of mild mass segregation for the faintest stars at the 
outskirts of the cluster. 

5.3. Degree of Stellar Stratification 

In order to parametrize the phenomenon of stellar stratifi- 
cation as it exhibits (or not) itself in the considered clusters, 
we assume as null hypothesis that all clusters are segregated. 
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TABLE 4 

Magnitude limits, the corresponding effective radii limits, 

AND THE DERIVED SLOPES Sstrat INDICATIVE OF THE DEGREE OF 
STRATIFICATION FOR SELECTED BRIGHTNESS RANGES FOR ALL THREE 
CLUSTERS OF OUR SAMPLE. 



Cluster 


mi4 


AmgM 


^cff 


Areff 






(mag) 


(mag) 


(arcmin) 


(arcmin) 


(mag/arcmin) 


NGC 1983 


12.0- 15.5 


3.5 


1.20-0.58 


-0.62 


-0.177 




15.5-22.0 


6.5 


0.58 - 0.94 


+0.36 


+0.055 




12.0-22.0 


10.0 


1.20-0.94 


-0.26 


-0.026 


NGC 2002 


11.0- 14.5 


3.5 


0.11-0.99 


+0.88 


+0.251 




15.0-23.0 


8.0 


0.49-1.23 


+0.74 


+0.093 




11.0-23.0 


12.0 


0.11 - 1.23 


+1.12 


+0.093 


NGC 2010 


13.0-14.0 


1.0 


1.24-0.10 


-1.14 


-1.140 




14.5-18.5 


4.0 


1.02-1.07 


+0.05 


+0.013 




19.0 - 23.5 


4.5 


0.89-1.51 


+0.62 


+0.138 




13.0-23.5 


10.5 


1.24-1.51 


+0.27 


+0.026 



Note. — Negative values of 5stiai or values close to zero indi cate non- 
segregated clusters at least within the specific magnitude limits (see j} l5.2) . The 
third row for every cluster refers to the whole observed magnitude range. 



but to a different degree. This is evident, as discussed above, 
from the r^ff vs. magnitude plots of Figure |7] This degree of 
stratification in each cluster may be, thus, expressed in terms 
of the brightness range of its stars, and the effective radius 
within which they are confined. As a consequence, in order to 
quantify stellar stratification as observed in every cluster, one 
may use the two primary output parameters of the effective 
radius method, i.e., the magnitude range and the correspond- 
ing effective radii of segregation. These parameters actually 
represent the ranges in magnitudes and radial distances that 
define the slope of the observed relation reff(msi4). Naturally, 
for a non-segregated cluster this slope, Sstrat, should be < 0, 
as it is derived from the graphs of Figure |2l Considering that 
a steeper positive slope of this relation represents a higher de- 
gree of stratification, we parametrize stellar stratification in 
the clusters by measuring 5strat = Areff/Am8i4 and quantify 
the phenomenon through this degree. 

As mentioned above, reff(mgi4) does not represent a system- 
atic relation throughout the whole observed magnitude range 
in every cluster, but, as seen in Figure |7J it shows differences 
depending on the selected magnitudes. For every cluster we 
select one indicative magnitude limit to separate the bright 
from the faint stars and measure the corresponding degree of 
stratification. This selection is based on the magnitude bin, 
where we observe dramatic changes in the relation reff(m8i4) 
in each cluster. As a consequence, we select two magnitude 
ranges, one spanning from the brightest magnitude bin down 
to the selected limit, and the other starting at this limit until the 
faintest magnitude bin, and we estimate their corresponding 
'Sstrat for all three clusters. These magnitude ranges, the corre- 
sponding Am8i4, and Argff, and the derived 5stiat are given in 
Table |4] along with the derived S'strat for the whole observed 
magnitude range for every cluster. The values of ^strat. given 
in this table, correspond to the degree of stratification for each 
cluster in every given brightness range. 

5.4. Results and Discussion 

From the ^snat measurements of Table |4] one can see that 
while indeed the bright stars of NGC 1983 are not segre- 
gated, stars with mgu ^ 15.5 mag show a definite trend 
in their relation reffCmgn), indicative of stratification with a 
degree of -1-0.06. On average, however, the cluster does not 
appear to demonstrate stellar stratification. We characterize. 



thus, NGC 1983 as partially segregated cluster. On the other 
hand, NGC 2002 shows a definite trend of larger reff for fainter 
stars within both selected magnitude ranges, as well as for 
the whole extend of observed magnitudes. It is a definite 
case of stellar stratification with a degree of stratification of 
-1-0.09. Finally, NGC 2010 is the most peculiar case, as the 
results for the few brightest stars are quite inconclusive, and 
for the fainter stars with 14.5 ^ ni%\4 ^ 18.5 mag the clus- 
ter appears slightly segregated, with a relation reff(m8i4) that is 
rather flat. On the other hand, the most prominent trend in this 
relation appears for even fainter stars with mg 14 18.5 mag 
and a degree of stratification S^tiAi — +0.14. However, if in- 
deed the cluster was segregated for this magnitude range, the 
observed trend should appear as a continuation of the (non 
stratified) relation reff(/«8i4) of the brighter stars. More specif- 
ically the Teff, where the trend starts for mg^ ~ 19 mag 0.'9) 
is smaller than that of the faintest non segregated stars of 
;7i8i4 ~ 18.5 mag with r^ff ^ 1^1. This behavior of the rela- 
tion reff(m8i4) coiTesponds to fluctuations in the spatial distri- 
bution of the faint stars around the center of the cluster rather 
than to true segregation. As a consequence, we characterize 
NGC 2010 as not segregated cluster. 

While the aforementioned results on three clusters cannot 
be used for a statistically complete interpretation of the phe- 
nomenon of stellar stratification in young LMC clusters in 
connection to their individual characteristics, it is worthwhile 
to search for any dependence of the degree of stratification 
to the structural and evolutionary parameters we derived for 
each cluster in our sample. Concerning the structural parame- 
ters shown in Table [3] all clusters are found with comparable 
concentration parameters, and tidal radii, r^, that are not far 
from each other. The most important difference we find in 
these parameters is the core radius, r^, of NGC 2010, the non 
(or less) segregated cluster, which is at least two times larger 
than those of the other two clusters. In addition, from the char- 
acteristics of Table |2] it is shown that this cluster is the most 
extended and least dense in the sample. Moreover, the CMD 
of the cluster (Figure|6l) shows indications of multi-age stellar 
populations, suggesting that NGC 2010 is a cluster with mul- 
tiple stellar generations. In such clusters a large fraction of the 
first generation of stars is lost early in the cluster evolution due 
to its expansion and stripping o f its outer layers result ing from 
early mass loss associated (e.g.. lD'Ercole et al.l2008h . In gen- 
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eral, the loose appearance of a low-density cluster with large 
rc, comes to the support of this speculation. While this can- 
not be verified with the present data, it would be reasonable to 
suggest that the lack of stellar stratification in this cluster will 
make the effect of early mass loss due to stellar evolution less 
destructive leading to longer lifetimes ( Vesperini et al. 2008*). 

As far as NGC 2002 is concerned, it is interesting to note 
that this cluster, being the most segregated in the sample, is 
actually the youngest (Table |2]l and located in a region, which 
appears to correspond to the emptiest general LMC field of 
the three (Table [3]). Naturally, these observations allow us to 
suggest that detection of stellar stratification may have been 
benefitted by the loose contaminating stellar field, and that 
the observed stellar stratification in NGC 2002 may be pri- 
mordial in nature due to the star formation process, or due to 
very early dynamical evolution. Indeed recent studies show 
that cool, fractal clusters can dynamically mass segregate on 
timescales comparabl e to their crossing-tim es, far shorter than 
usually expected (e.g.. lAllison et al.l2009bl) . The exact nature, 
though, of the observed segregation cannot be fully under- 
stood without detailed kinematic information about the clus- 
ter members, or deeper photometry that will provide the com- 
plete IMF of the cluster Finally, the observed partial stel- 
lar stratification in NGC 1983 is driven by the fact that the 
brightest stars of the cluster are distributed in a wide range of 
distances from the cluster center In general, the behavior of 
Teff as function of stellar magnitude cannot be connected to 
any specific difference of its characteristics from those of the 
other clusters, apart from the fact that the cluster belongs to 
the richest region in field populations, and it is a rather ellip- 
tical cluster 

6. CONCLUSIONS 

We present a coherent comparative investigation of stellar 
stratification in three young LMC clusters with the application 
of a robust method for the assessment of this phenomenon on 
deep HST imaging. In a cluster where stellar stratification 
occurs, the segregated brighter stars are expected to be more 
centrally concentrated than the non-segregated fainter stars. 
Stellar stratification can, thus, be investigated from the de- 
pendence of the radial extent of stars in specific magnitude 
(mass) ranges on the corresponding r nean magnitude (mass) . 
According to this notion, in Paper I jGouliermis et al.ll2009h 
we developed and verified the efficiency of a robust method 
for the assessment of stellar stratification in star clusters. This 
method is based on the calculation of the mean-square radius, 
the effective radius r^ff, of stars in different magnitude ranges, 
and the investigation of its dependence on magnitude as indi- 
cation of stellar stratification in the cluster. 

In the present study we apply the effective radius method 
for the detection and quantification of stellar stratification in 
young star clusters in the LMC. We select three clusters ob- 
served with the high resolving efficiency of HST/ ACS, specif- 
ically NGC 1983, NGC 2002 and NGC 2010, on the basis 
of the differences in their appearance, structure, stellar con- 
tent, and surrounding stellar field. Our photometry delivered 
complete stellar catalo gs do wn to 111^4 ~ 22 - 23 mag for all 
three observed fields (§ 12. lb . The stellar surface density maps 
of the observed regions, constructed from star counts on the 
photometric catalogs, demonstrate that all three clusters are 
morphologically quite different from each other. NGC 1983 
appears to be a centrally concentrated, rather elliptical com- 
pact cluster, NGC 2002 a compact, spherical stellar concen- 
tration, and NGC 2010 a large, loose and rather amorphous 



cluster (§[TB. 

The applic ation of bot h EPF (lElson. Fall & Freeman 1 1987h 
and King's (iKingI [T96 2') models on the stellar surface den- 
sity profiles of the clusters allowed us the measurement of 
the background stellar field density of the regions, and the 
estimation of the core radii, r^., of the clusters. While these 
clusters do not seem to be tidally truncated, a tidal radius, 
is estimated, indicative of the limiting radius of each cluster 
(§ El). The latter is essential for the successful application 
of a random field subtraction technique for the decontamina- 
tion of the stellar samples of the clusters from their surround- 
ing young stellar ambient and the general LMC field (§ 14.31 ). 
We derive the ages of the clusters from the stellar populations 
comprised within the of the clusters and we find ages of 
r ~ 28 Myr, 18 Myr and 159 Myr for NGC 1983, NGC 2002, 
and NGC 2010 respectively (§1421). 

We finally apply the effective radius method for assess- 
ing stellar stratification in the clusters and for its quantitative 
study (§ |5]l. We bin the stars according to their magnitudes 
in the F814W filter, and we estimate the corresponding effec- 
tive radii, reft, for every cluster in our sample. We then plot 
the derived radii versus the corresponding mean magnitude 
and we investigate the functional relation, reff(m8i4), between 
these parameters. With our method it is shown that stellar 
stratification behaves differently in every cluster NGC 1983 
appears to be partially segregated, since its brightest stellar 
content does not appear to be centrally concentrated, while 
its fainter stars show a mild dependency of reft towards lower 
values for fainter magnitudes. The results on the bright stars 
of NGC 2010 show evidence of lack of stratification, since 
the brightest stars are located quite far away from its center. 
As far as the faint stars are concerned they show a rather flat 
reff(OT8i4) relation, with no evidence of segregation in this re- 
lation for the faintest stars. Consequently, this cluster is char- 
acterized as not segregated. Finally, NGC 2002 is found to 
be well segregated for both its brighter and fainter stars. It is 
a clear case of proof of stellar stratification with the effective 
radius method. 

We propose the slope 5strat = Areff/Am8i4, measured for 
stars in selected magnitude ranges, as well as in the whole 
magnitude range, as the most appropriate parameter for the 
quantification of the degree of stratification, and we use it to 
characterize the phenomenon as i t is o bserved in the clusters. 
From the derived values of 5stiat (§ 15. 3b we find that NGC 2002 
is the most strongly segregated cluster in the sample with a 
degree of stratification 5stiat — +0.09, NGC 1983 is a partially 
segregated cluster with Sstm — +0.06 for the faintest stars with 
m^i4 ^ 15.5 mag, and NGC 2010 is not segregated with in- 
dicative S'stiat — +0.03 (§ 15.4b . Finally, it should be noted that 
the extension of this study to a larger sample of LMC clus- 
ters observed with HST will provide a more complete picture 
of the phenomenon of stellar stratification. Naturally, the de- 
velopment of such a comparative scheme, which will include 
star clusters in the whole extent of the evolutionary sequence, 
requires a comprehensive set of ACS and, in the near future, 
WFC3 observations. 
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